function dfout = rhs(t,u,L)

% Nx = length(u);
% uhat = fft(u);
% kap = (2*pi/L)*[-Nx/2:Nx/2-1];
% kap = fftshift(kap');
% duhat = i*kap.*uhat;
% du = ifft(duhat);
% dfout = -du;

Nx = length(u);
uhat = fft(u);
kap = (2*pi/L)*[-Nx/2:Nx/2-1];
kap = fftshift(kap');
duhat = i*kap.*uhat;
dduhat = -(kap.^2).*uhat;
du = ifft(duhat);
ddu = ifft(dduhat);
dfout = -u.*du + .01*ddu;